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Abstract 

The mixing boundary-value problem for infinite one-dimensional chain of harmonic 
oscillators on the half-line is considered. The large time asymptotic behavior of solutions 
is obtained. The initial data of the system are supposed to be a random function which 
has some mixing properties. We study the distribution fit of the random solution at time 
moments t € R. The main result is the convergence of fit to a Gaussian probability 
measure as t —>• 00 . The mixing properties of the limit measures are studied. 
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1 Introduction 

We consider the following infinite system of harmonic oscillators on the half-line: 

u(x, t) = (Al — m 2 )u(x, t), ieN, t > 0 , ( 1 . 1 ) 

with the boundary condition as x = 0 : 

it(0, t) = F(u(0, t)) — m 2 u(0, t) — 7 ?i( 0 , t) + w(l, t) — w(0, t), t> 0, (1.2) 

and with the initial condition (as t — 0 ) 

u(x,0) — u 0 (x), u(x, 0) = v 0 (x), x > 0. (1.3) 

Here u(x,t) Gl 1 , m > 0, 7 > 0, A i denotes the second derivative on Z 1 : 

A L u(x) = u(x + 1) — 2 u(x) + u(x — 1), x G Z 1 . 

If 7 = 0, then fjl.ip - fll.2p is a Hamiltonian system with the Hamiltonian functional 

H(w, u) := - ^2 (|h(:r )| 2 + | u(x + 1) — u(x )| 2 + m 2 \u{x)\ 2 ^j + P(u( 0)), (1.4) 

Z a :>0 

where, by dehnition, P(q) — f F(q)dq, q G M 1 , stands for the potential energy of the 
external force. To prove the existence of solutions to problem fjl.ip ~ fjl.3p we assume that 
P = 0 or 

P G C 2 (R), P(q) ~> Too as 1^1 —» oo, (1.5) 

so P(q) > P 0 f° r ah q with some Pq G M. 

Write Y(t) = (F°(t), F 1 ^)) = (u(-,t),ii(-,t)), Y 0 = (l^do 1 ) = (mo(-)Wo(-)) • We assume 
that the initial state Fo(^) belongs to the Hilbert space T-i a ,+ , ol G M, consisting of real 
sequences, see Dehnition 12.11 below. The existence and uniqueness of the solutions Y(t) is 
proved in Appendix A. 

To prove the main result we assume that F(q) = —Kq with k > 0. Moreover, on the 
coefficients m, k, 7 of the system we impose some restrictions. Namely, if 7 > 0 we assume that 
7Ti > 0 or m — 0 and If 7 = 0, it is assumed that kG (0, 2). The initial state Y 0 (x) 

is supposed to be a random element of the space P n ,+ , ot < — 3/2 , with the distribution /i 0 . 
We assume that fi 0 is a probability measure of mean zero satisfying conditions S2-S4 below. 
In particular, we assume that the initial measure /i 0 satisfies a mixing condition. Roughly 
speaking, it means that Y 0 (x) and Y 0 (y) are asymptotically independent as \x — y\ —> 00 . 

For a given f G R, denote by fit the probability measure on T-i a ,-7 giving the random 
solution Y(t) = ()) to problem (11 . ID (II.31) . Our main objective is to prove the 
weak convergence of the measures /i t on the space 7-L a ,+ with a < —3/2 to a limit measure 
Ho o , which is an equilibrium Gaussian measure on T-i a ,+, 

fjit —r Hoo as t —> 00 . ( 1 . 6 ) 

This means the convergence of the integrals j f(Y) n t (dY) —> j f(Y)n oa (dY) as t —> 00 
for any bounded continuous functional / on . 
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For one-dimensional chains of harmonic oscillators in the whole line, the convergence to 
equilibrium have been established by Boldrighini et al. in [T] and by Spolin and Lebowitz in 
[23j. Ergodic properties of one-dimensional chains of anharmonic oscillators coupled to heat 
baths were studied by Jaksic, Pillet and others (see, e.g., mm)- 

In [21 S], we studied the convergence to equilibrium for the systems described by partial 
differential equations. Later on, similar results were obtained in [3] for d -dimensional harmonic 
crystals with d > 1, and in [6] for a scalar field coupled to a harmonic crystal. The convergence 
(11.61) was derived for the harmonic crystals in the whole space 7L d , d > 1, and in the half-space 
Z := {x = (xi,... ,Xd ) € Z d : x\ > 0} with zero boundary condition (see [5J and [7], resp.). 
In the present paper we prove the convergence (11.61) only in the one-dimensional case but with 
the boundary condition of the form (11.21) . 

We outline the strategy of the proof. Using the technique from 0E], we derive the con¬ 
vergence (11.61) from the assertions I and II: 

I. The family of measures fi t , t > 0, is weakly compact in T-i a ,+ with a < —3/2 . 

II. The characteristic functionals converge to a Gaussian functional, 

AtO) : = f exp(i(y, T) + ) fi t (dY) -> exp{-iQ 00 (T, T)}, t ->■ oo, (1.7) 

Here T = (T°, T 1 ) e S := S®S , where S denotes a space of real rapidly decreasing sequences, 
<y,tf> + = E E Y\x)^ i \x) for Y = (y 0 ,^ 1 ) e U a + and T = (T°, T 1 ) €5, Q 00 ('F, T) is 

i=0,l x>0 

a quadratic form. 

To prove I we derive the uniform bound (15.ip below for the mean local energy with respect 
to the measure fit , t > 0, and apply the Prokhorov Theorem. 

To prove II we study the asymptotic behavior of the solution Y ( t ) and obtain (see 
Lemma [6.71) that 


(Y{t),*)+ ~ (C/o(f)y 0 ,n*), t oo, (in mean), (1.8) 

where Uo(t) is a solving operator to the mixing problem (I2.2j) — (I2.4D with zero boundary con¬ 
dition, the vector-functions are expressed by T e S (see formula (16. 19j) below). 

To prove (ll.8[) we decompose the solution u(x, t ) into two terms: u(x, t ) = z(x, t ) +q(x, t ). 
Here z(x,t) is a solution of (II. ip with zero boundary condition and satisfying the initial 
condition (11.31) . i.e., z(-,t)) = Uo(t)Yo, q(x,t) is a solution of (11. ID with zero initial 

condition and some boundary condition (see (12.71) below). The existence and behavior of the 
solution z(x, t ) were studied in [7]. In Section [2731 we state the basic results on z(x, t ) (see also 
Appendix C). Therefore, the next step is to study the long time behavior of q(x,t ), x > 0. 
In the case when x > 1, we obtain the following formula for q(x, t ), 

t 

q(.z> t) — J K(x,t — s)q(0, s) ds, x € N, (1.9) 

0 

where the kernel K(x,t) is defined in (13.51) below. The properties of this kernel are studied in 
Appendix A. In particular, we prove that K(x,t ) satisfies the following bound 

^^(1 + x 2 ) a \K(x, t)\ 2 ~ (7(1 + t)~ 3 as t —> oo, for any a < — 3/2. (1-10) 
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Therefore, the main step in the proof of (jl. 8 j) is to study the behavior of q(x,t ) with x = 0. 
By (11.21) and (II. 9p . q(0,t) evolves according to a Volterra integro-differential equation of a 
form 


t 

q(0,t) = F(q(Q,t)) - (1 + m 2 )q(0,t) —'yq(0,t) + J K(l,t — s)q(Q, s) ds + h(t), t > 0, ( 1 . 11 ) 

o 

with h(t) = z(l,t). To study the long time behavior for the solutions of Eqn (II. lip , we put 
F(q) = —Kq with k > 0, apply the Fourier-Laplace transform t —> co to the solutions q( 0, t) 
and study the analytic properties of < 7 ( 0 , 0 ;) for complex values co G C (see Appendix B). 
These properties allow us to obtain the following bound for the solutions q( 0, t) of Eqn (II. lip 
with h{t) = 0 : 

\q(0,t)\ < C(l + t)“ 3/2 , t> 0 . 

Applying this estimate together with (11.91) and ( 11 . 101 ) . we obtain the asymptotics for q(x,t) 
in mean (see formula (16.161) below) 

(<?(•, t),ip)+ ~ (U 0 (t)Y 0 , K°) + , t ->• 00 , 

where ij> £ S , the vector valued function K° is defined in (16.151) below. This implies the 
asymptotics (II. 8 p which plays the crucial role in our convergence analysis for the statistical 
solutions of the problem (ll.ip — (jl.3D . 

The dynamics of the equations with delay has been extensively investigated by many au¬ 
thors under some restrictions on the kernel K(l,t) . For details, we refer to the monograph 
by Gripenberg, Londen and Staffans [D]. The stability properties of linear Volterra integro- 
differential equations can be found in the papers by Murakami [T9], Nino and Murakami [2D] , 
Kordonis and Philos [16]. Note that in the literature frequently the long time asymptotical 
behavior of solutions is studied assuming that the kernel has the exponential decay or is of one 
sign. However, in our case, by ( 11 . 101 ) . the decay of K(l,t) is like (1 + f )~ 3 / 2 . 

In recent years the nonlinear equation of a form ( 11 . 111 ) with a stationary Gaussian process 
h(t) and with a smooth (confining or periodic) potential P(q) = — f F(q)dq has been in¬ 
vestigated also extensively, for example, the ergodic properties of such equations were studied 
by Jaksic and Pillet in [12] . and the qualitative properties of solutions were established by 
Ottobre and Pavliotis in [21] , In the present paper, we prove the convergence to equilibrium 
for a linear model. However, we do not assume that the initial distribution of the system is a 
Gibbs measure or absolutely continuous with respect to a Gibbs measure. Therefore, the force 
h(t) = z(l,t) in Eqn (II. lip is non-Gaussian, in general. 

The paper is organized as follows. In Section [2] we impose the conditions on the model 
and on the initial measures fi 0 and state the main results. The limit behavior for solutions of 
Eqn (II.lip is studied in Sections [3] and [4] The compactness of the measures family {fi t , t G M} 
is proved in Section O In Section £0 we establish asymptotics ( 1 1 . 8 p and complete the proof of the 
convergence (11.61) . In Appendix A we study the properties of K{x,t ) and prove the existence 
of the solutions. The behavior of q(0,t) is studied in Appendix B. Appendix C contains the 
results on the solutions z(x,t ) . 
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2 Main Results 


2.1 Phase space 

We assume that the initial date Y 0 belongs to the phase space H a ,+ , aft, defined below. 

Definition 2.1 (i) £^ + = £ 2 + (Z+), aeR, is the Hilbert space of sequences u(x), x > 0 , 
with norm 

IMI 2 ,+ = H^l 2 ^) 20 < °°> ( x ) : = (! + kl 2 ) 1/2 - 

x>0 

(ii) H a ,+ = £ 2 < 8 ) £ 2 is the Hilbert space of pairs Y = (■ u , v ) of sequences equipped with 

norm 

\\Y\\l,+ = IM| 2 ,+ + IM| 2 ,+ < oo. 

H. a ^ is equipped with the Borel a -algebra B(fH a ,+) ■ 

(Hi) £ 2 a = ^(Z 1 ) is the Hilbert space of sequences with norm ||u || 2 = Y \u(x)\ 2 (x) 2a < oo . 

x£7j 

In particular, £\ = £ 2 . Write H a £ 2 a ® £ 2 a , a 6 l. 

Theorem 2.2 Let 7 , m > 0 and condition \1.5\) hold, and let Y 0 G H a . + , a G R . Then the 
mixing problem 11.11 1- 0731) has a unique solution Y(t) G (7(M, H ai+ ) ■ The operator U(t ) : 
Y 0 —>■ Y(t) is continuous on H a ,+ ■ Moreover, for any T > 0 , sup| t | <T \\U(t)Y 0 \\ aj+ < B, 
where the constant B depends on ||lo||a,+ and T. 

Theorem 12.21 is proved in Appendix A. The proof is based on the following representation 
for the solution u(x,t ) of the problem (ll.ljMjl.3p : 


u(x,t) — z(x,t) + q(x,t), x > 0, ( 2 . 1 ) 

where z(x, t ) is a solution of the mixing problem with zero boundary condition, 

z(x, t) = (A l — m 2 )z(x, t), igN, t > 0, (2.2) 

z(0,t) = 0, t > 0, (2.3) 

z(x,0) — u 0 (x), z(x, 0)—v o (x), igN. (2.4) 

Therefore, q(x,t ) is a solution of the mixing problem with zero initial condition, 

q[x, t) = (A l — m 2 )q(x, 1), iGN, t > 0, (2.5) 

9(0, t) = F(q( 0, 1)) - m 2 g(0, t ) - jq( 0, t ) + q( 1, t ) - g(0, t ) + z( 1, t), t > 0, (2.6) 

<7(x, 0) = 0, g(x,0) = 0, igN, (2.7) 

9 ( 0 , 0 ) = u 0 ( 0 ), 9 ( 0 , 0 ) = w 0 ( 0 ). ( 2 . 8 ) 


In Section [231 we state the results concerning the solutions z(x,t ). In Sections [3] and 51 we 
study the behavior of the solutions q(x, t) . 

To prove the main result we assume that F(q) = —nq , with k > 0. Moreover, on the 
nonnegative constants m, 7 , k of the system we impose the following condition C. 

C. (i) If 7 yl 0 , then either m ^0 or m — 0 and k 7 ^ 0 ; 

(ii) if 7 = 0, then k G (0, 2). 
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2.2 Random initial data 

We assume that the initial date Y 0 (x) = (l^°(x), Yf(x)) is a measurable random function with 
values in (fH a ,+, B(fH a ,+)) • Denote by /i o a Borel probability measure on TL a) + giving the 
distribution of Y 0 . The expectation with respect to /U 0 is denoted by E. 

We impose the following conditions S1-S4 on . 

51. The measure fi 0 has zero expectation value, E(Y 0 (x)) = f Y 0 (x)/u 0 (dY 0 ) = 0, 

52. The measure /i 0 has hnite variance, supE(|l'o(x)| 2 ) < eo < oo. 

x>0 

Write u 0 = /i 0 {lo G TLa^ : ho(0) = 0} . The expectation with respect to is denoted by 
E 0 . On the measure we impose conditions S3 and S4. 

53. The correlation functions of the measure u 0 are denoted by 


Q 0 0,0 = E 0 (Dg( x)Yq (^O)) x , x ' > 0) i,j = 0,1. 

In particular, the matrix Q 0 (x,x') = (Qq (x, x'))ij =0j i vanishes if x = 0 or x' = 0 . Moreover, 
for every x G Z 1 , 

lim Ql\x + y,y) = q 0 \ x ). 

y—H-oo 

Here qf (x) are correlation functions of some translation invariant measure v with zero mean 
value on the space TL a . By definition, a measure v is said to be translation invariant if 
u(T h B) = v(B) for B G B(TL a ), h G Z 1 , where T h Y 0 (x) = Y 0 (x — h ) for igZ 1 . 

To formulate the last condition on z/ 0 , denote by A an open interval of and by cr(A) 
the a -algebra on TL a ,+ generated by Y 0 (x) with x G A. Define the Ibragimov mixing 
coefficient of a probability measure v on TL a ,+ by the rule 


V>{r) = 


sup 

A,B 

dist(M, B) > r 


sup 

A G a(A), B G a(B) 
u(B) > 0 


u(AHB) - v(A)v(B)\ 

US) 


(2.9) 


Definition 2.3 A measure v is said to satisfy the strong uniform Ibragimov mixing condition 
if ip(r) —> 0 as r —> oo . 

S4. The measure uq satisfies the strong uniform Ibragimov mixing condition with mixing 
coefficient ip , and 

+oo 

J ip l ^ 2 [r) dr < oo. (2.10) 

o 

The examples of the initial measures z/o satisfying S3 and S4 are given in [7]. 

Lemma 2.4 Let conditions S1-S4 hold. Then qf G i 1 . Moreover, for any <f>, T G "Ho.+ , the 
following bound holds 


\(Qo(z, z'), *(z) ® * (/)>+| < C||$|| 0 , + ||tf ||o,+. 


( 2 . 11 ) 
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Proof Indeed, by [[TO:, Lemma 17.2.3], conditions SI, S2 and S4 imply 

\Qo ( z , z ')\ < Ceoip 1/2 {\z~£\), z,z'eN, i,j = 0,1. (2.12) 

Hence, condition (12. 10p implies (12.13(1 . 

y l/2 (\z\) < C\ < oo for all z6N, (2-13) 

z'eN zez 1 

and 

Y \Qo( z i Ol < Ci < oo for all zJ £ N. (2.14) 

zGN 

Here the constant C\ does not depend on z, z' £ N. Estimates (12.131) and (12.141) and the Shur 
lemma imply the bound (12.111) . Moreover, by condition S3 and the bound (12.121) . 

\q l 0 J (z)\ <Ce 0 g> 1/2 (\z\), zeN. 

Hence, q £ i x . This implies, in particular, that qf £ (^(T 1 ) . Here and below q{9) = 
F x ^.g[q(x)\ denotes the Fourier transform w.r.t. x £ Z 1 , 

q(6) = Y ), 6* £ T 1 = M72t rZ 1 . ■ 

x EZ 1 

For a probability measure v on denote by 0 the characteristic functional (Fourier 

transform) 

z>(T) = J exp (i(Y, T) + ) v(dY), T £ 5. 

A measure v is called Gaussian (with zero expectation) if its characteristic functional has the 
form b(T) = exp{ —|Q('L, T)} , £ S , where Q is a real nonnegative quadratic form in S . 


2.3 The mixing problem with zero boundary condition. Convergence 
to equilibrium for problem ( 12.211 —( f2~4ll 

In this section we study the problem (I2.2D - 02.4D . The results of this subsection were proved in 
[7j. To state these results, write Z(t) = Z(x,t) = (z(x,t), z(x,t)). 


Lemma 2.5 (see Lemma 2.1 in J7f) Assume that a £ M . Then the following assertions hold. 

(i) For any Yq £ TL a ,+ , there exists a unique solution Z{t) £ C(M.,'H a ,+) to the mixed problem 

(ii) For any t £ M ; the operator t/ 0 (t) : Y 0 K y Z{t) is continuous on TL a ,+, and 


sup \\Z(-,t)\\ 2 ^ 

\t\<T 


<C(T)\\Y 0 


The proof is based on the following formula for the solution Z(x, t) of problem (12.2D (12.41) : 


z*{x,t) = y Y^+^ XjX ^ Y o^ x '^ x G z +’ c- 15 ) 
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where Z(x,t ) = (Z°(x,t), Z 1 (x,t)) = (z(x,t), z(x,t)), the Green function Q t + (x,x') is 

1 


Q tt+ (x,x') := G t {x-x')-G t (x + x'), Gt(x) = -£-l e lxe Q t {9)d9, x G Z 1 , (2.16) 

Z7T Jj 1 


with 


cos 0(0)£ 


sin 


ftW = (C?W)U.= 1 C0 * s l V 


, 0(0) = v/2-2cos0 + m 2 . (2.17) 


In particular, 0(0) = 2| sin(0/2)| if m, — 0 . We see that Z(0, t) = 0 Vt, since G t (—x) = ^(x) 
(see formulas (12.161) and (]2.17jl ). 

Definition 2.6 We define u t as the Borel probability measure on TL a ,+, which gives the dis¬ 
tribution of solution Z(t ) to the problem A2.2\) i.e., 

vfiB) = u 0 (Uo(-t)B), where B G B(fH a ,+), tel. 

Lemma 2.7 The correlation functions of measures u t converge to a limit, 


Ql J (x,x')= / Z l (x)Z :, (x') vfidZ) —* Q^x, x'), t —> oo, x,x' eZ|. 


(2.18) 


The correlation matrix Qoo(x,x') = (Q^x, x'))jj =0 has the form 

Qoo(x,x') = qoo(x - x') - goo(x + x') - qoo(-x - x') + qoo(-x + x'), x,x'el} + . (2.19) 

Here qoo(x) = (x) + (x), x G Z 1 , where in the Fourier transform, we have 

qW) = \m) + c{e)UO)c T m = \ 


1 ( q°o° + $V“ 2 (0) € - 


Qoo(d) = |sgn(0)(C'(0)go(0) -q 0 (9)C T (9)) 

. f % 10 -C 0 2 (^o°° + 0o n 

- IWV) ^ -0 2 (0)g o oo-^i 0 2 (0)(g o lo -g o 01 ) 


'Qo 1 + 0o° 0 2 M° + 0o n " 


( 2 . 20 ) 


idere 0 G T 1 if m ^ 0 and 0 G T 1 \ {0} if m = 0, (-) T denotes a matrix transposition, the 
functions q'f , i, j — 0,1, are introduced in condition S3, 

C(0) = ^ ^ , 0(0) = V2 — 2 cos 9 + m 2 . (2.21) 

Denote by Q(( 0 ('h, \h) a real quadratic form on S defined by 

= E E ( 2 . 22 ) 

*,l=o,i x,x'ezi 
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Remark 2.8 Given T = (4°, T 1 ) G £^ + x i \ + , introduce an odd sequence \k 0 (x ), 16 Z 1 , 
such that 'ko(x) = 'k(x) for x > 0, \k o (0) = 0 and \k 0 (x) = — \k(— x) for x < 0. Note that 
Q 00 (x,x') is odd w.r.t. x and i'gZ 1 . Moreover, 

i 

S£,(*,*) = <to(* - A ® 9'W)) = E E 9 ~( 1 (2-23) 

*,j=0 a^a/GZ 1 

In the case m/0, gT G , by Lemma 12.41 and formulas (12.200 . Then, Young’s inequality 
yields 

|< 2 ^(vk,vk)| <G||g 00 ||,i||^ 0 ||^<G 1 ||g 00 ||,i||vk||^ + , 

i.e., Q£o(\k, T) is continuous in "Ho ,+ = *g, + x £l, + • Here as before || • ||o (|| • ||o,+) stands for 
the norm in £ 2 = £ 2 or in TLq (in £q + and in Lfo,+ , respectively). 

In the case m = 0, (f)(6) = 2|sin(0/2)|. Note that \4> 0 (8)\ < C'(^)|sin0| for if; e S. 
Denote t\ /4>& = {if) G ^ 0 ,+ : e Y. 2 (T 1 )} with norm ||V , ||i/^,+ : = IK 1 + 1 / < / , ( 0 ))V’o( 0 )IIl 2 ( T i) ■ 

Therefore, for any T = (T°, T 1 ) G 7Yi/^, + := + x £ 2 0)+ , 

Q^(T, T) < IkolUi {C.WF-^l/J)) * 4r«||g + C 2 \\*o\\l) < C 3 \\^\\l^ + + G 4 ||T 1 || 2 + , 

by (12.201) and (12.231) . Here and below by F~ l we denote the inverse Fourier transform. 

Theorem 2.9 Let conditions SI, S3, S4 hold and let a < —1/2 if m ^ 0 and a < — 1 if 
m — 0 . Then the following assertions are fulfilled, (i) The measures u t = Uo(t)*vo weakly 
converge as t —* 00 : 

u t > z/qo as t ^ 00 on £he space % ai+ . 

/n) The limit measure is Gaussian, with zero mean value and with correlation matrix Q 00 
defined in l\2.19\) . 

(Hi) The uniform bound holds, 

supE 0 ||C/ 0 (f)Yo||a + < 00 . (2.24) 

t> 0 

Moreover, lmp^E 0 |([/o(t)Y 0 , ^) + | 2 = for an V ^ G S . 

(iv) The limit measure ^ satisfies a mixing condition w.r.t. the group U 0 (t), i.e., for any 

f,geL 2 (n a , + ,u «,), Inn J f(U Q (t)Y)g(Y)u 00 (dY) = f f(Y)v 00 (dY) J g(Y)v 00 (dY) . 

If m/0, the results of Lemmas 12.71 and Theorem 12.91 follow directly from [7]. In the case 
m — 0, the proof needs in some modification. In Appendix C we prove the bound (12.241) for 
m > 0 which implies the compactness of the measures family {u t ,t G M} , and outline the 
proof of Theorem 12.91 

2.4 Main theorem 

Definition 2.10 (i) We define fit as the Borel probability measure on TL a ,+ given by the rule 
fit(B) = Ho(U(—t)B ), where B G B((H a ,+) and t G M. 

(ii) The correlation functions of the measure fi t are defined by 

Q)((x,x') = E (Y\x,t)Y j (x',t)y i,j = 0,1, x,x’ G N. 
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Here Y l {x,t ) are the components of the random solution Y(t) = (V°(-, t), t)) to problem 


Denote by T) a real quadratic form in S of the form 

Qoo(d','h) = Q«(n®,n«), ( 2 . 25 ) 

where Qoo is defined by (12.221) . 11^ from (16.191) . 

Remark 2.11 For any 16 5, 11^ G FU/cf,^ ■ This follows from l\6.19\) and Remark Id. ,51 
Therefore, by Remark\2fE and formula \2.25\) . Q 0 o(d / , T) is bounded for all T G FLo,+ ■ 

The main result of the paper is the following theorem. 

Theorem 2.12 Let conditions S1-S4 on the initial measure /i 0 hold and the constants m, 7 , n 
satisfy condition C. Then the following assertions are fulfilled. 

(i) The measures p, t weakly converge as t —> 00 to a limit measure on the space Li a)+ 
with a < —3/2 . Moreover, the limit measure is Gaussian with zero mean and the correlation 
matrix Qoo . 

(ii) The correlation functions of p t converge to a limit, and for T G <S, 

lim E| (Y(t), v h) + | 2 = Q 00 ( v h, T). 

t—¥ OO 

3 Large time asymptotic behavior for q(x, t) 

In Sections [3] and |U we investigate the behavior of the solutions q(x,t ) to the problem (j2 ,5j) — 
(12.8D as t —> 00 . At first, we study the properties of q(x, t ) with x > 1 in the Fourier-Laplace 
transform. I 11 next section, we study the asymptotics of q(0,t). 

Definition 3.1 Let \q(t)\ < Ce Bt 
formula 

q(u) 

The Gronwall inequality and conditions on F(q) imply standard a priori estimates for the 
solutions q(x, t) , x > 1. I 11 particular, there exist constants A, B < 00 such that 

t)\ 2 + \q(x, t)\ 2 ) < Ce Bt as t —> 00 . 

ccEN 

This bound is proved in Appendix A (see formula (17.191) 1. Hence the Laplace-Fourier transform 
with respect to t -variable, q(x,t ) —> q(x,ui) , exists at least for Imca > B and satisfies the 
following equation 

(— A l + m 2 — uj 2 )q(x, oj) = 0, igN, Imw > B, (3.2) 

by (12.71) . Now we construct the solution of (13.21) . Note first that the Fourier transform of the 
lattice operator —A l is the operator of multiplication by f>o(9) = 2 — 2 cos 9, and f>o(9) G 
[0,4] . Second, the following lemma holds (see Lemma 2.1 in [15]). 

Denote T m := [—\/4 + m 2 , —m\ U [m, \/4 + m 2 ]. 


The Laplace-Fourier transform of q(t ) is given by the 

+OO 

J e lut q(t) dt , Imw > B. (3.1) 

0 
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Lemma 3.2 For given uj G C m = C \ T m ; the equation 

2 — 2 cos# = oj 2 — m 2 (3.3) 

has the unique solution (we denote it by 9(oj) ) in the domain {6 G C : Irn 6 > 0, — 7r < 
Re# < 7r}. Moreover, 0(oj) is an analytic function in C m . 

Since we seek the solution q(-,t) G f 2 a + with some a , q(x,u> ) has a form 

q{x,u) — q{D,u)e l6 ^ x , x > 0. 

Applying the inverse Laplace-Fourier transform with respect to a;-variable, we write the so¬ 
lution q(x, t) in the form 

+00+ifl 

q(x, t) — — [ e~ lult q(0,uj)e ie( ' UJ ^ x dco, i6 N, t > 0, fi > a, (3.4) 

27r J 

— OO+Z/i 


and q(-,t) G (( + . Moreover, the left-hand side of (13.4)1 does not depend on /i. The integral 
in (13.41) is understood in the sense of principal value. Write K(x,cu) = e*°A) x an d put 

+oo 

K(x,t) — — [ e~ lU}i K(x, co) div, igN, t > 0, /i > a. (3.5) 

27T J 

— OO+2/X 


Then (13.41) implies the following representation 

(q(x, t), q(x, t)) — I K(x, t — s)(q(0, s), g(0, s)) ds, iGN. (3.6) 

Jo 

To study the large time behavior of q(x,t ), we will move down the contour of integration in 
(13.4)1 and in (13.5)1 . In Appendix A, we study the analytic properties K(x,cu) for w G C and 
obtain the following result. 

Theorem 3.3 For every iGN, K(x,t ) = 0 for t < 0 . Moreover, for any a < —3/2, the 
following bound holds, 


J2{%) 2a \K{x,t)\ 2 <C(l + t)~ 3 for t> 0. 




In particular, 

\K(l,t)\ < C(l + f)" 3/2 , t> 0. 
Corollary 3.4 For any a < —3/2, Id. 61) and \3. 71) imply 


(3.7) 


(3.8) 


1/2 


XA) 2 “(W I ’*)| 2 + W 1 ' 4 )! 2 ) ) < C l(l+t-s) 3/2 (k(0, s)| + |g(0, s)|) da, t> 0. 


xGN 
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4 Large time asymptotic behavior of g(0, t) 

Let F(q) = —nq with k > 0. Write q(t) = q(0,t). By (13.61) . Eqn (I2.6[) becomes a linear 
Volterra integro-differential equation with the kernel K(l,t ), 

t 

q(t) — — (k + 1 + m 2 )q(t) - yg(f) + J K(l,t - s)q(s) ds + z(l,t). (4.1) 

o 

At first, we study the solutions of the corresponding homogeneous equation (with z(l,t) = 0 ), 

t 

q(t) — —(k + 1 + m 2 )q(t) - ^q(t) + J K(l,t - s)q(s) ds, (4.2) 

o 


with the initial data 

9(0) = «o(0) = 9o, 9(0) = u o (0) = Pq. (4.3) 


Remark 4.1 For m > 0 , set 


'7m ■ — 


+ \/4 + m 2 


m 


(4.4) 


Note that 7 m = max \4>'(0) \, where (f)(6) = \/2 — 2 cos 6 + m 2 . Moreover, by direct calcula- 

06[—7T,7r] 

tion, we obtain 

r+o o 

/ K{l,s)ds = K(l,ti) = e*W = 'r i m . 

Jo 

In particular, if m ^ 0 ; G (0,1) . If m = 0, then = 1. 

Applying the Fourier-Laplace transform (13.11) to the solutions q(t) of Eqn (14.2)) . we obtain 

D(uj)q(oj) = —icoqo + go7 + Po for Irn uj > B, 


where, by definition, 

D(u) :=—u 2 + k + 1 + m 2 — icj'y — K(l,u), K(l,u) = e 10 ^. (4.5) 

Write N(oj) := [i)(a;)] _1 . Then q(u>) = N(uj ) (—iu>q 0 + g 0 7 + Po) • The analytic properties of 
D(uf) and N(uf) for oj e C are studied in Appendix B. In particular, N(oj) is analytic in the 
upper half-space, i.e., with Irn u > 0. Denote 


m 


i 

27T 


-\-oo-\-ifi 

J e~ iut N(u) du ;, 


— 00 + 2/4 


t > 0, 


with some p > 0. 


(4.6) 


The following proposition follows from Theorem 18.51 (see Appendix B). 
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Proposition 4.2 Let constants m, 7, k be nonnegative and satisfy condition C. Then 


\N^{t)\<C(l+t)~V\ j = 0,1,2. (4.7) 

Moreover, the solution q(t ) of the problem has a form 

q(t) = N(t)(p 0 + 7 g 0 ) + N(t)q 0 . 

Hence, the following bound holds, \q{t)\ + \q(t)\ < C( 1 + |t|)~ 3 / 2 (|<?o| + bol) f or an U t>0. 

Remark Let 7 = 0 . If k — 2 or m ^ 0 and k — 0, then |7V^(t)| < C( 1 + t ) -1 / 2 . For 
m — k — 0 with any 7 > 0 , N(t) = (7 + l )” 1 + o(l) as t —> 00 (see Theorem 18.51 below). 

Corollary 4.3 Let condition C hold. Denote by S(t ) a solving operator of the Cauchy problem 
<ra- Then the variation constants formula gives the following representation for the 
solution of problem fl 4-l\ ), (O-- 


q(t) 

q(t) 


= s(t ) 


9o 

To 


+ /V) 

Jo 


0 

z(M-r) 


dr, t > 0. 


Evidently, S'(O) = /. The matrix S(t) , t > 0 , is called the resolvent or principal matrix 

solution for E,„ 0. Moreover, the matnx S{t) has a form ( £[*> /(/ 

Proposition \f.S\ implies that | S' (7) | < C(1 + t)~ 3 / 2 , and for the solutions of \f.l\ ) the following 
bound holds: 


\q{t)\ + \q(t)\ < Ci(l + t) 3/2 (|g 0 1 + |po|) + C 2 / (1 + s) 3/2 \z(l,t - s)| ds for t > 0. (4.8) 

Jo 


5 Compactness of measures family 

The compactness of the measures family {p, t , t > 0} follows from the Prokhorov Theorem and 
the bound EH] (see below), since the embedding TL a ,+ C TLp,+ is compact if a > /3. 


Lemma 5.1 Let a < —3/2 and conditions SI S3 and C hold. Then 

supE||r(-,t)|£ )+ < C < oo. 


t> o 


(5.1) 


Proof By the representation (12.ip . the solution Y(t) has a form Y(x,t ) = Z(x,t ) +X(x,t), 
where X(x,t) = (q(x,t), q(x,t)). Hence by the bound (12.241) . to prove the bound (15.11) it 
suffices to verify that 

supE||X(-,t)|| 2 < C < oo. (5.2) 


t> o 


Applying formula (13.61) and bound (13.71) . we have 

E||V(V)II«, + = E|.Y(0,f)| 2 + E5>> 


2 a 


iSN 
2 


K(x, t — s)X (0, s) ds 


< C sup E|X(0, r)\. 

re[0,i] 


(5.3) 
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Using the estimate (14.8p . we obtain 

E|X(0,t)| 2 < E(|g 0 | 2 + M 2 ) + CsupE|z(l, r)| 2 , t> 0. 

r> 0 

On the other hand, sup rgR E|,s(l, r)| 2 = sup TgR Q°°(l, 1) < C < oo by the bound (12.24}) . 
Hence, 

supE|X(0, t)\ 2 < C < oo. (5.4) 

T> o 

Estimates (15.31) and (15.41) imply the bound (15.21) . since 

supE||U(-,t)|| 2 < supE||Z(-,t)|| 2 +su P E||X(-,t)|| 2 < C < oo. ■ 

t> 0 t> 0 t> 0 


6 Asymptotic behavior of Y(t) in mean 


6.1 The estimates of q(0,t) in mean 

Set q(°)(x,t) = q(x,t ), q^(x,t) = q(x,t ) , To formulate the asymptotic behavior for 

q^(0,t) in mean, introduce the following notations. By definition, (see (12.16P ) 

Gi(y,t):= (e,“ + ( z,v),Silt(z,yj) = (Qf(z-y)-g?(z + y),Q',\z-y)-g?(z + y)), (6.1) 

y, z £ Z 1 , i — 0,1, leR. In particular, G*(0, i) = 0 for any z and t. Let G fiy) denote 
the vector valued functions defined as 

Too Too 

G%)= j N( s )G{(y,-s)ds= J N^(s)G°fy,- S )d S , y e Z 1 , j=0,l, (6.2) 

0 0 

where the function iV(s) is introduced in (14.61) . G{(y, s) is dehned in (16.11) . 

Proposition 6.1 Let conditions SI S3 and C hold, Y 0 £ 'H a ,+ ? an d 9(0, t) be the solution 
to problem flUT Then 

g W (°T) = (U 0 (t)Y 0 ,G J ) + + r j (t), j = 0,1, where E\rj(t)\ 2 < C(1 + t) -1 , (6.3) 

the vector valued functions G J , j =0,1, are defined in Ii6. &) . 


Proof First, Corollary 14.31 implies that 


E 


9(0,t) 
9(0, t) 


/ sw U."-n 


2 


dr 


< C(l+t)~ 3 


Second, S(r) 


0 

Z(l,t~ T ) 


N(r) 

N(t) 


z(l, t — r). Moreover, 
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Further, the bound (12.241) gives 


E(z(l,t - 7i)z(l,t - r 2 )) 


< C supE 0 |z(l, s)| 2 < C < oo. 

s£ R 


Hence, the bound (14.71) yields 


E 


+ OO 


S(r) 


0 

Z(l,t~ T ) 


dr 


+oo 


< C'supE|^(l,s)| 2 [ / (1 + r) 3 / 2 dr ] <Ci(l + £) 1 . 


This implies the representation (16.31) with j — 0, since by (16.11) . 

^(1, t — r) = (Uo(t)Y 0 (-), G?(-, —t)) + . 

The bound (16.31) with j = 1 is proved by the similar way. ■ 

Remark 6.2 (i) By (j2.16j) and (I2.17P . G l z (y,t) is odd w.r.t. y E Z 1 . Moreover, G*(t/, i) = 0 
for z — 0. Also, the Parseval identity gives 


2tt 


7T 


G I *(-,0IIS = ^ / |GUd,t)| 2 dd = - / (|£f(d)| 2 +|^ 12 


sin 2 (^d)d6 l , Vz G N. (6.4) 


Hence, if m / 0, then ||G*(-, t)\\l < C < 00 uniformly in z E Z 1 and t E R. If rn — 0, 
||'Gl (-,t)||g < C < 00 uniformly on z and t, and 

l|G°(.,0llg < C, + c 2 r S1 " 2( ;( ) do < c, + CM- 

Jn sill (9/2) 

In particular, G i(9,t) = 2i sin 6^ cos((p(9)t), sm(tp(9)t) / <p(9) j , 6 E T 1 . Then formulas (12.171) 
and (16.41) imply 


7r 



(cos 2 ((j)(6)t) + 


sin 2 (</>(6*)t) 

W) 



—7V 

(ii) G J (-) is odd. Moreover, by the bounds (14.71) and (16.51) . G j E Ho , since 


(6.5) 


+ OO 


+ OO 


\G j (-) 


lo < 


|iV ( - ?) (s)|||G t j ) (-, —s)|| 0 ds < C / \N®(s)\ds<oo. 


( 6 . 6 ) 


In Fourier transform, 


+00 

G j (9) — 2 i sin 9 J N^\s) (cos(0(#)s), — s'm((fi(9)s)/4>(9)Jds, 9 E T 1 . 

0 

Therefore, G J E Hy# : = ly^ x i 2 , where E 2 /^ = {i/> E l 2 : ^/(j) E L 2 ( T 1 )}. Note that 

= ^ if ™^0. 
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Denote by U' 0 (t) the operator adjoint to U 0 (t): 

(Y, = {U 0 (t)Y, *>+, for Y G n a ,+, $65, t G R. (6.7) 

In other words, 

1 

for *6 5, te R, yez\. 

i=0 a;>0 

In particular, (see (16.1(1 ) G°(z/,l) = ((7o(f)Yo)(z/) with Yo(x) = (5iz,0), where <5i x denotes the 
Kronekker symbol. 

Corollary 6.3 /i) For G J defined in H6.2 1) . we hare 

supE|(I/o(f)ho, G J ) + | 2 = supE|(y 0 ,Co(t)G J ) + | 2 < C < oo. (6.8) 

ZgR tgR 

Indeed, since U^fiG^y, r) = G°(j/,r + t), then 

sup ||(7g(t)G : '||o < C < oo, (6.9) 

tgR 

by 16. 61) . Hence, H6.S 1) follows from the bounds \2.11\) and 16.9 1) . 

/ii) The representation 16.1) . Lemma\2/^ and the bound \6. 5)1 gire the following convergence, 

E(g«(0,t)g w (0,t)) -)• Q^(G*,G J ) as t ^ oo, (6.10) 

where the quadratic form is defined by \2.22fi . The r.h.s. of Id. 101) is defined by Re- 
marks UT7H and 1 6. 21 (ii). 

6.2 The large time behavior of q(x,t ), x G N, in mean 

Let K7(x, z/) j = 0,1, stand for the vector-valued functions, 

+00 +00 +00 

K J (x,z/) = J /l(x,s)^I/o(-s)G^(z/)ds = J j K(x, s)N b) {r)G\{y,-s - r) ds dr, (fill) 
0 0 0 

i 6N, i/EZ 1 , K(x,s ) is defined in (13.5(1 . G J is introduced in (16.21) . Note that K ? (x, y) is 
odd w.r.t. y , K J (-, y) G %<*,+ with a < —3/2, and K7 (x, •) G T£o,+ for any x . 

The formula (13.61) and Proposition 16.11 imply the following lemma. 

Lemma 6.4 (i) The solution q(x,t) , x > 1, of problem \2.tfi - \2.8\) admits the following 
representation 

q {3 \x,t) = (U 0 (t)Y 0 ,K J (x, ■))+ + rfix,t), j — 0,1, t > 0, (6.12) 

where E||rj(-,f)||o + < C(1 + t)^ 1 as t —> oo, K7 is introduced in Id. 111) . There ||r||o + = 

E,e*K*)l 2 - 

(ii) The correlation functions converge to a limit, 

lim E (q( l \x,t)q ( j\x , ,t) \ = QLi K*(x, TKdx'-)) for any x,x' G Zi, i, j = 0,1. 

t->oo V / 
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Proof At first, by (13.61) and (16.31) . 


q^(x,t) — / K(x, t — s)q ( ' 3 \0, s) ds — / K(x,t — s)(Uo(s)Y 0 , G 3 (-)) + ds + r'Jx,t), (6.13) 


where x G N, E||r'(-, t)||g + < C(1 + t) 1 . Indeed, by (16.31) and (13.71) . 

n\r'(;t)\\l + = E [ K{-,t-s)rj(s)ds q+ < (J \\K (-,t - s)|| 0 ,+ y / IE|r i (s)| 2 ds^j 


< c( / (1 + t — s) 3 ^ 2 (l + s) ds) <Ci(l + t) 




Second, the hrst term in the r.h.s. of (16.131) has a form 


K{x, t - s)(U 0 (s)Y 0 ,&(■))+ ds = ( U 0 (t)Y 0 , K*(x, •))+ + r"(x, t ), 


where, by (16.111) . 


r»+oo 


r"(-,t)= / K(x,s)(U 0 (t - s)Y' 0 , G J (-)) + ds. 


Since G J G fg,+ x ^o,+ j the bounds (I2.24H . (13.71) and (16.81) yield 
+ 00 

E||rJ(-, i)HS,+ < (J ||if(-,s)lk + (Eo|{C/o(i-s)ro.Gn.)) + | 2 ) 1/2 d S ) 2 <C(l + «r I . (6.14) 


Hence, (16.131) and (16.141) imply (16.121) with rj(x,t) = r){x,t ) +r"(x,t) . Finally, the item (ii) 
of Lemma 16.41 follows from the representation (16.121) and Lemma 12.71 ■ 

Remark If we set K( 0,cu) := e l0 ^ x \ x=o = 1, then K(0,t) = 5 0t and K4(0,y) = G J (y). 
Hence, we put K4(0, y) = G J (y) , y G Z 1 . Then the representation (16.31) follows from (16.121) . 

For any ip e S, denote 

/•+oo 

Kj(y) = (K J '(-,S/),^(-))+= / (K(;s),iP) + (U'(-s)&)(y)ds, ye Z 1 . (6.15) 

Jo 


Remark 6.5 IVofe first that G J (0) = 0 (see K6.2 j) L Hence K J (x, 0) = 0 /or any x G N . 
Then K^(0) = 0 . Also, is odd. Moreover, K 3 ^ G f 2 ^ x f 2 (see Remark WM (ii)). 

Corollary 6.6 (i) For any ip G S', j = 0,1 7 

(q {:> \-it),ip) + = (U 0 {t)Y 0 ,K 3 ll) ) + + r j (t), t > 0, where E|r,-(f)| 2 < C(1 + f) _1 . (6.16) 

(ii) For any ip G S , 

supE| (Uo(t)Y 0 , K/) + | 2 = supE|(F 0 ,f/'(t)K;) + | 2 < C < oo. (6.17) 

te R te R 

(in) For any ip,x 6 S, E ((#(-,f),^) + (g 0) (-, t),x)+) -> 2^o( K ^ K i) as t ^ oo . 
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Proof The item (i) follows from (16.121) . Now we verify the item (ii). 

By definition (16.151) . s), i)j) + U' Q {t — s)G^{y) ds. Hence, (16.91) implies 

the following uniform bound, 

f+OO 

sup||f/o(t)K^|| 0 < C(^) / ||A'(-,s)|| 0 .+ ds < cx). (6.18) 

te R Jo 

Therefore, (16.171) follows from the bounds (12. lip and (16.180 . Finally, (16.161) and Lemma 12.71 
imply the item (iii) of Corollary 16.61 ■ 

6.3 Proof of Theorem 12.12 

Write W(x, y) = ej5 xy + K J (x, y) with eo = (1, 0), e.\ — (0,1), j = 0,1, x, y G 7} + . Given 
T = (T°, T 1 ) G S , define the vector valued functions n$(?/), y G , 

l l 

n*(») := »),**)+ = t(a) + 5 ](K J (-,!,).>py)) + = Mv) + Y 2 K Uv) (6.19) 

3=0 3=0 

with K lp(y) from (I6.15p . The formula (12.11) and Lemma [6.41 imply the following lemma. 

Lemma 6.7 For the solution Y(t ) = U(t)Y" 0 = (u(x,t),u(x,t)) of the problem QHy ( COj) , 
the following asymptotics holds, 

(U(t)Y o y(x) = (U 0 (t)Y 0 ,W(x,-)) + + r j (x,t), j =0,1, (6.20) 

where rj(x,t) from and E||r,-(-,£)||o )+ < 67(l + t) _1 . Hence, for any T = (T°,T 1 ) e S, 

< Y(t ), T) + = (U 0 (t)Y 0 , n^) + + r(t), (6.21) 

where r(t) — X^)=o( r i('^)> ^)+ > am d E|r(t)| 2 < 67(1 + t)^ 1 . 

We return to the proof of Theorem 12.121 The compactness of the measures family {nt,t G 
M} is proved in Section [5j To end the proof of item (i) of Theorem 12.121 it suffices to check 
the convergence (11.71) of characteristic functionals for /i t ■ By the triangle inequality and the 
equality (12.251) . 


Ee i<r(t),*)+ _ e -i<2oo(*,*) 

< 

w ( JV _ p hf/ohho.n*> + A 



V ) 


+ 


Eoe^ t/ °^ y °’ n ’H + 


_ g-fQSHmp.n*) 


( 6 . 22 ) 


Applying (16.211) . we estimate the first term in the r.h.s. of (16.221) by 

1/2 


E 


(Y(t),V)+-(Uo{t)Y 0 ,Il s ,) + <E|r(t)| < (E|r(t)| 2 j < 67(1 + t)“ 1/2 -)• 0 as t -y oo. 


The convergence of the characteristic functionals h t (n^) = E 0 (exp{i(C/ 0 (^)lo, n^) + }) to a 
limit as t —> oo follows from Theorem 12.91 and Remark 12.111 This proves the convergence 
(11.71) . Corollary 16.61 and Lemma [2.71 imply the item (ii) of Theorem 12.121 
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7 Appendix A: The behavior of K(x,t) as t —> oo 

We first study the analytic properties of K(x, oo) := e ld ^ x (where Q{uS) is introduced in 
Lemma [3. 2 1) applying the technique of Komech, Kopylova and Kunze |Toj . developed by them 
for discrete Schroclinger and Klein-Gordon equations in the whole space Z 1 . 


7.1 Properties of e l0 ^ x for w G C and x 6 N 


As before we denote T m := [—\/4 + m 2 , — m] U [m, \/4 + m 2 ] . For x G N, C, the function 
K(x,oo ) = e l0 ^ x has the following properties I III. 

I. Let oo G C m — C\ r m . Then Im0(w) >0. In this case, K(x,oo) exponentially decays in 
x . Hence K(x,oo) is an analytic function in the complex oo -plane with the values in the class 
+ = + • The behavior of K(x,oo) for |cu| —» oo is following. 

Let Imcu — pt — const > 0. Then lm.6(oo) — > +cx) and Re#(u;) — > ±7r as Rew —* ±oo. 
(If Imcu = fjL < 0, then Im0(w) —> +oo and Re$(uz) —> =F7r as Recu —» ±oo .) Moreover, for 
any K > 0 there exists a positive constant C = C(K) < oo such that 

| e »0(i/+i/d I < f or any i/gl: \v\ > \/4 + m 2 , and any p G M : Ml < K. (7.1) 

1 1 v 2 

Indeed, put z = 1 — ( oo 2 — m 2 )/ 2. Then, by (13.31) . 

_ cos ^ ^ sin 0 = 2 + A 2 - 1 = 


s/z 1 - 1 ’ 


(7.2) 


where \/— 1 is a complex root and its branch is chosen by the condition Im 0(oo) > 0 . 
In particular, this condition implies that sgn(Im \J z 2 — 1) = — sgn(Im z) = sgn(Imo; 2 ) for 
Imu 2 7 ^ 0. Moreover, \Jz 2 — 1 > 0 for w G (—oo, —\/4 + m 2 ) U (\/4 + m 2 , +oo). Let us 
verify (17.11) for oo = v G M: |i/| > \/4 + m 2 . For such values oo , Re(#(z/)) = ± 7 r and 
\+ 2 — 1 = +(2 — v 2 + m 2 ) 2 /4 — 1 > 0. Hence, by (17.21) . 


& ie(u) _ g—im o{y) — 


2 — (z / 2 — m 2 ) — a/(2 — z / 2 + m 2 ) 2 — 4 


as zd —>■ oo. 


(7.3) 


In the general case, the bound (17.11) can be proved similarly. Moreover, by direct calculation, 
we have |e* e ^)| < C|a ;| _2 as |u;| —* oo. 

Lemma 7.1 The bound 17.11) implies that there exist constants C, B < oo such that 

'^^(x) 2a \K(x,t)\ 2 < Ce Bt for any a Gl, t > 0. (7.4) 

igN 

Proof Let us apply the inverse Fourier-Laplace transform. Then for some // > 0 . any x G 11 
and t > 0, 

+oo +ifi -foo 


K(x, t ) = 


2t r 


e -wt e i0( w )* = 


0 flt 

2n 


e -ivt e iO(v+in)x dv ^ 


— OO+2/l 


where oo = u + ip. Note that 


ji8(y+ip)x 


__ g—Im 6(y+in) 


= : < 1 for any izGl, and /r > 0. 
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Hence, 




< 


xGN 


n 2 fit 

47r 2 


+ 00 


(Y,Y) x dv ■ 


xGN 


Therefore, (17.41) follows from (17.11) . since 


£<*>*“ (O' s 

xsn 


C 


, ^ oo, for a < 0, 


Ewkj < cx;w<c 1T ^ 


tcEN 


zGN 


1 - Qv* 


C 


\v\ 


, v —y oo, for a > 0. 


II. Let w 6 1 and |u;| G (m, \/m 2 + 4) . Then cn 2 — m 2 G (0,4) and Im0(a;) = 0. Put 
9(ou ± iO) = lim £ _ ) . + o 9(ou ± ie) . For every ieN, the following pointwise limit exists 

K(x, oj ± ie) —> K(x, oj ± iO), £ —> +0. 

Moreover, 1 9{oj ± ie)\ < C for any e . Hence, for any a < — 1/2 , 

|| K(x,oj ± *0) — K(x,lu ± 7£)|| 2 )+ = ^^(x) 2a \K(x,cu ± iO) — K(x,oj ± ie)\ 2 —» 0 as e —» +0, 


:cEN 


by the Lebesgue dominated theorem. 

Note that 0(u;) = —9(ou) for oj G C m = C \ T m . Hence, 


K(x,oj — iO) = K(x,oj + iO) for a; G (—\/4 + m 2 , —m) U (m, V4 + m 2 ), x G N. 

III. Now we discuss the behavior of K(x,oj) near the points ±m, ±\/4 + m 2 . Eqn (13.31) 
implies 


e *6»(o;) _ cos q^ _|_ i s in #(a;) = 1 — -(a; 2 — m 2 ) + iy oj 2 — m 2 — -(w 2 — m 2 ) 2 , a; G C m . 
The Taylor expansion implies 

e *0H = 1 + i\/a; 2 — m 2 — i(a; 2 — m 2 ) — Yoj 2 — m 2 ) 3//2 + ..., a; —* ±m + iO, (7.5) 

2 8 

where sgn(Re \Joj 2 — m 2 ) = sgn(Rea;) for Imw > 0 . This choice of the branch of the complex 
root \/ a; 2 — m 2 follows from the condition \m9[oj) > 0. Hence, for iGN, 


+oo 


e i0(u) x = 1 + J2(u 2 - m 2 ) j / 2 R j (x), lu ->• ±m + iO. 
i=i 


(7.6) 


Here polynomials W(x) are of a form /P (x) = c ,i xk w ihh some coefficients c(gC, j G N. 

k=\ 

For example, R l {x) = ix , R 2 (x) = —x 2 /2 , R 3 (x) = —i(4x 3 — x)/24. 

The behavior of K(x,cu) near points oj = ±\/m 2 + 4 is similar. Namely, 

e ie ( w) = -1 + i(m 2 + 4 - u; 2 ) 1/2 + ^(m 2 + 4 - w 2 ) - ^(m 2 + 4 - (n 2 ) 3/2 + ..., (7.7) 

2 8 
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where \Jm 2 + 4 — u 2 is the complex root and we choose the branch of this root such that 
sgn(Re \Jm 2 + 4 — to 2 ) = sgn(Reiw) by the condition linden) > 0. Hence, for x G N, 


e i6 ^ )x = (H) I (l-ttv / m 2 + 4-aj 2 -(m 2 + 4-a;V/2 + .. 
If m — 0, then the Taylor expansion gives 


, , ,2 , ,4 

^ CU . CO UJ 

= 1 - +iuj(l -+ .. 

° v 8 128 


as u> —* ±Vm 2 + 4 + 10. (7.8) 


ca —)■ 0, (7.9) 


and e 10 ^ = —1 + iy/ A — a; 2 + ... as a; —* ±2 + 10 . Therefore, 

£,i8(u)x _ i -)_ _ uj 2 x 2 /2 — la; 3 (4x 3 — x)/24 + ... as uj —> 0 

e i8(uj)x _ ^ 1)^(1 — ixy /4 — uj 2 — (4 _ uj 2 )x 2 /2 + ...) as a; —» ±2 + 10 

The representations (17.611 . (17.8(1 and (17.101) imply the following lemma. 
Lemma 7.2 (cf JT5l Lemma 3.2]) Let m — 0 . Then for every ARM, 

N 

e i8(uj)x UJ J R J (x) + Tn(u},x), uj —> 0 , 

1=1 


(7.10) 


where ||rjv(a?, -)IU,+ = C ) (|a;| Ar+1 ) f or a < —3/2 — At. 

Lei m/0. Then for every N e N, 

N 

e i8{u) x = x + ^(w 2 — m 2 y/ 2 R j (x) + rjv(w, x), a; —>• ±m + 10, (7.11) 

l=i 

where ||rjv(w, -)||a,+ = CK(a; 2 — m 2 )^ JV+1)//2 ) for a < —3/2 — At. Moreover, 

N d r 

D r u ( e «(a»)*j = — (^ - m 2 ) J '/ 2 iT» + f N {u,x), uj^±m + 10, (7.12) 

l=i 

where ||rjv(a;, •)|| a ,+ = 0({uj 2 — m 2 )( JV + 1 )/ 2 - r ) f or a < —3/2 — At. The similar representation 
holds for lu —> ±y/m 2 + 4 + 10 . 

Indeed, the bound (17.111) follows from the following representation for remainder (see for¬ 
mula (17.61) ) 

JV+l 

r N (co, x) = (uj 2 — m 2 y N+1 ^ 2 ^ b k (uj)x k , 

k=l 

where b k (uj) are uniformly bounded for uj —> ±m. In particular, 

£ i8(uj)x _ x + y/u 2 — 7 n 2 R 0 (uj , x) as uj —> ±m, 
where sup ||i2o(o;, -)l|a,+ < C < oo for a < —3/2. 
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7.2 Proof of Theorem 13.3 


By the properties I III, K(x,t ) = 0 for t < 0. To prove (13.7|) . we use the properties I III 
and apply the technique from m- At first, we rewrite Ii(x,t ) in the form 


K(x, t) = — 
v ; 2 tt 


e zut K(x,u)du = -/ e lult K(x, u) du, x E N, 

27t J 
r 


Im uj=fi>0 

where T = {|u;| = R : R > \/4 + m 2 } and T is oriented anticlockwise. Since K(x,u) is 
analytic in C\T m , we can change the contour of integration into T m s , where the contour T m>£ 
surrounds of T m and belongs to the £-neighborhood of T m (r m , £ is oriented anticlockwise). 
Taking e —y 0 , we obtain 


K(x,t) = — I e 


—iut 


K(x, u + iO) — K(x, oj — zO)) du 


= - I e 

7T „ 

f 


— iidt 


Im K(x, u + iO) du = — I e lult P±(x, u) du. 


± 3 =1 


Here Pf(x, u) := Q (co>)Im/t (x,<n + iO) , where £■ (u), j — 1,2, are smooth functions such 

that )T) C 7 ± ( w ) = 1, u; G R, supp C (9(±m), suppC^ C 0(±\/4 + m 2 ) (Here O(a) 
±d 

denotes a neighborhood of a point u = a). If rn — 0, we introduce Cj (Pi) instead of ^ 
(P± , resp.) with suppCi C 0(0) . By the property III, 

||Pi ± (-,o;)|| Qi+ = 0(\u =Fm| 1/2 ) if m ^ 0, ||Pi(-,w)|| a)+ = 0(|w|) if m = 0, 

HTf(-,^)lk+ = 0(|a; =F \/4 + m 2 |) 

for any a < —3/2. Therefore, using the lemma of Jensen and Kato | [T4I Lemma 10.2] or the 
Vainberg lemma [251 Lemma 2], we obtain 


e luJt P±(x,u) du = 0(t 3 / 2 ), t —> oo, j = 1,2, for any a < — 3/2. (7-13) 


The bound (13.7(1 is proved. 


7.3 Existence of solutions 

Lemma 7.3 Let a E R, m, 7 > 0, and let P satisfy the condition U.5\) . Then the following 
assertions hold, (i) For every Yq E TL a .+ , the mixing problem M. JD -/ TO) has a unique solution 
Y(t ) € C(M, TL a ,+) ■ Moreover, the operator U(t ) : To 1 —)■ Y'(t), tel, is continuous on FL a ,+ , 
and there exist constants C, B < 00 such that 

||T(t)|| a>+ < C'e s *||r 0 |U+ fortER. (7.14) 

(ii) For Y 0 E P 0 ,+ , 

H(Y(t)) +7 [ |u(0,s)| 2 ds = H(Y 0 ), tel, (7.15) 


where H(Y(t )) is defined in ( l.f). In particular, if 7 = 0 , the energy H(Y(t )) is conserved 
and finite. 
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Proof of Lemma 17.31 We use the representation (12.11) . u(x,t ) = z(x,t ) + q(x,t). By 
Lemma [2751 and formula (13. 61) . to prove the existence of u(x,t ) , it suffices to prove the ex¬ 
istence of the solutions q(t) = q(0,t) to the problem (14.2p ~ fj4.3j) . Indeed, first, rewrite (14.21) in 
the equivalent integral form 

t S t S 

q(t) — I I 7F(T,q(r)) drj ds + J ^J z(1,t) dr — yg(s)j ds + C 0 + Cfi, t > 0. (7.16) 

0 0 0 0 

with 7F(t, q) := -(1 + m 2 )q + F(q) + f*K(l,t- s)q(s ) ds , C 0 = q(0 ) = q 0 , C x = q{ 0) = p 0 ■ 
By the contraction mapping principle, the solution q(t) of Eqn (17.161) is unique on a certain 
interval t G [0, e) with some e > 0 depending on the initial data qo,Po ■ Hence, by (13.61) . 
the solution q(x, t ) of problem (I2.5D - 02.8D is unique on a certain interval t G [0, e) with some 
e > 0 depending on Y 0 . The existence of z(x,t ) is stated in Appendix C. This implies the 
existence of the local solution u(x,t ) = z(x,t ) + q(x,t ) for sufficiently small t. This local 
solution can be extended to the global solution using the a priori estimate (17.141) . Now we 
verify (17.141) . Indeed, by (14.21) . 


- (m\ 2 + (m 2 + l)\q(t)\ 2 )+P( g (t)) + 7 / \q(s)\ 2 ds 


(bol 2 + {m 2 + 1)M 2 ) + P{qo) + f q(s)(z(l,s)~ j K(l, s - T)q(r)drj ds. (7.17) 


o 


0 


De&ie M(t ) = sup (|q , (s)| 2 + |q , (s)| 2 ). Then the equality (17.171) and the bound (13.81) for K(l,t) 

o <s<t 

yield 

t t 

M(t) < Ci + J \/M(s)\z(l, s)\ds + C 2 J M(s ) ds, t > 0. 

o o 

Applying the Gronwall-Bellman integral type inequality [22] we have 

t 

M(t) < e C2t (-\fC~i + ^ f \z(l,s)\e C2S,2 ds) 2 , t > 0. 


Since |z(l,t)| < C( 1 + |i|) JV ||^o||a,+ ( see Appendix C), we obtain the a priori bound 

\q(t)\ + \q(t)\ < Ce m \\Y 0 \\ a , + 


(7.18) 


with some constants C, B < oo . By Lemma [7. II and formula (13.61) . 

[^2(x) 2a (\q(x,t)\ 2 + |g(x,t)| 2 )) 1 < C'ie Bt |[r 0 ||a,+- (7.19) 

zGN 


Since ||P(t)|| 2 i+ < ||g(-,t)||q, + + ||gM)||! I+ + ll-(-^)l| 2 ,+ + P(-^)l| 2 ,+, the a priori bound 
(17.14j) follows from (17.191) and (19.51) . ■ 

Remark 7.4 Let F(q ) = —nq with k > 0, ho G 7-^o,-t- • Then the energy H(Y(t )) is nonneg¬ 
ative and finite, H(Y(t )) < H(Y 0 ) by 17.151) . 
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8 Appendix B: The analytic properties of D(uS) for wGC 

In this section, we study the properties of D(oo) and N(oo) = for wfC. 

Denote by C + (C_ ) the upper (lower, resp.) half-plane, C± = {u 6 C : ±Imw > 0} . 

Lemma 8.1 (i) N(oo) is meromorphic for to G C m — C \ Y m . 

(ii) For any K > 0 , there exists a constant C(K) > 0 such that for any p G K. such that 
\fi\ < K , we have 

\D(y + ifj)\ > C{K)v 2 , v —7 oo. (8.1) 

(Hi) (D(a ;)) _1 is analytic for oo G C + . 

Proof (i) It follows from the formula (14.51) and from the property I of K(x,oo) (see Ap¬ 
pendix A) that D{oj) is analytic for oo G C m . 

(ii) For any (3 > 0, there exists a constant Kg > \/4 + m 2 such that \D(oo)\ > C\oo\ 2 for 
any oo G C with |o;| > Kg. Hence, (D(cj )) _1 exists for large |u;| . Moreover, item (ii) follows 
from (14.5p and (17.ip . Hence, for any /iGR such that |/i| < K there exists C\ = C\ (K) > 0 
such that \N(i; + ip) \ < C\v~ 2 as v —> oo . 

(iii) Assume that D{ooq) = 0 for some ooq G C + . Introduce a function q*{t ), tGl, such 

that q*(t) = for t > 0, and q*{t) vanishes for t < 0. Then q*{t) is a solution of 

Eqn (14.2p with the initial data (1, — iooo) ■ Denote by Y*(x) := (uo(x),vo(x)) the initial data 
for the problem (ll.lD - (ll.2p such that Y*(x) = 0 for x G N and F*(0) = (1, — ico 0 ). Then the 
solution of the system (I2.2p - (l2.3p with the initial data F* vanishes, i.e., 

[/ 0 (t)F* = (z(-,t), z(-,t)) = 0 for any t > 0 and iGN. 

In particular, z(l,t) = 0. Hence, g*(t) is the solution of (14.ip for t > 0. Then q*(x,oo) := 
e * 6 »(w)x ~^(g, a; ) _ ^ e ie{u})x j^ _ Therefore, the solution of the problem (II . ID — (II.2f) with the 

initial data F* is of the form u*(x,t) = q*(x,t) = e * 0 (wo)® e -iwot , ^ e pj ; t > 0. On the other 
hand, the Hamiltonian H{Y(t )) (see (11.41) ) becomes 

77(u*(-, f), u*(-, t)) = e 2tlmuj °H(Y*) for any t > 0, where H(Y *) > 0. 

Since Im ooq > 0 and F* G "Ho,+ , this exponential growth contradicts the energy estimate 
(17.151) . Hence, D(oo) 0 for any oo G C + . ■ 

Corollary 8.2 7/7 = 0, then D(oo) ^ 0 for any oo G C_ . 

Indeed, if 7 = 0, then D(x,oo) = D(x,oo ), because 6{oo) = — 6{oo) for oo G C m . Then the 
assertion of the corollary follows from item (iii) of Lemma 18.11 

Now we study the invertibility of D(oo + iO) for w G 1. 

Lemma 8.3 Let oo G R, and the constants 7 ,m,K satisfy condition C. Then D(co-\-iti) 7 ^ 0 
for any 00 G M . 

Proof Step T. Let 00 G R and |(u| > \/4 + m 2 . Then Re ^(a; + f0) = ± 7 r. Therefore, 

D(u + iO) = - 00 2 + k + 1 + m 2 - ioo'f + e - Im0 ^ +i °), with Im 6{oo + iO) > 0. 
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Hence, Im D(uj + *0) 7 ^ 0 iff 7 7 ^ 0. On the other hand, 

Re77(cu + *0) = — cu 2 + k + 1 + m 2 + e - Im 0 O+ l °) < K — 2 for |cu| > \/4 + m 2 , 

and Re7)(cu + i0) = n —2 for a; = ±\/4 + m 2 . Moreover, Re7)(cu + iO) —>■ —00 as |cu| —* cx) . 
Hence, for |cu| > \/4 + m 2 , Re7)(cu + iO) 7 ^ 0 iff k < 2. Therefore, for such values u u, 

D(oj + *0) 7 ^ 0 iff 7 7 ^ 0 or 7 = 0 and n < 2. If 7 = 0 and k > 2 , then there exist two 

points ±cuo (cuo > \/4 + m 2 ) such that 77(±cu 0 + 10) = 0 . 

Step 2: Let m 7 ^ 0 and cu G (—m, 0) U (0,m). For such values cu, Red(cu + *0) = 0 and 
e id(u+t 0 ) e (e l 0 (o), e * 0 (±"*+*o)) _ ( 7 ^^ 1 ) with 7 m defined in (14.41) . Hence, 

Re D(u> + iO) = —cu 2 + n + 1 + m 2 — > « for |cu| < m, 

and Re 7)(cu + iO) = k if cu — ±m . Therefore, D(oj + iO) 7 ^ 0 for any |cu| < m , since k > 0 . 

Step 3: Let cu G (—\/4 + m 2 , —m)U(m, \/4 + m 2 ). Then Im#(cu + i0) = 0 and Re0(cu + iO) G 

(—7T, 0) U (0, 7 r). Moreover, sign (sin 0(cu)) = sign a;. Hence, for such values cu and for m 7 ^ 0 , 

Im D(oj + *0) = —CU 7 — sin0(cu) = — sgn(cu) (jcu |7 + \/cu 2 — m 2 \J\ — (cu 2 — m 2 )/4j. 

If rri = 0, 77 (cu + *0) = k — cu 2 /2 — iu ^7 + y/l — cu 2 /4^ . Therefore, for such values cu, 

Im 7?(cu + iO) 7 ^ 0 for any k, 7 > 0 . 

Step 4- For cu = ±\/4 + m 2 , e t0 ( w+ ®°) = —1, and 77 (cu + iO) = k — 2 7 = rp\/4 + m 2 . Hence, 

77(±\/4 + m 2 + 10) 7 ^ 0 iff 7 7 ^ 0 or 7 = 0 and k 7 ^ 2. 

For cu = ±m , e* 6 ^" 1 "* 0 - = 1, and 77 (cu + iO) = k =f R/m . Hence, for m/0, 77 (±m + iO) 7 ^ 0 

iff 7 7 ^ 0 or 7 = 0 and k / 0. For m = 0, 77 ( 0 + iO) = k 7 ^ 0 iff k / 0. Lemma 18.31 is 

proved. ■ 

Corollary 8.4 Let rn,^,K > 0. If constants m, 7 , k do not satisfy condition C, then there 
exist points oj G M in which 77 (cu + iO) vanishes. Indeed, if m — k — 0 , ihen 77(0 +10) = 0 
for any 7 > 0. If 7 = k — 0, i/ien 77(±m + iO) = 0 for any m > 0 . 7/7 = 0 and k — 2, 
then D(±\/m 2 + 4 + *0) = 0 for any m. 7/7 = 0 and k > 2, then there exists a point 
lo 0 > \/4 + m 2 snch that D(±.ujq + iO) = 0 . 

Now we study the behavior of D (cu) in the neighborhood of the special points cu = ±m 
and cu = ±\/4 + m 2 . 

For cu —» ±\/4 + m 2 + ?'(), the representation (17.71) holds. Hence, in the neighborhood of 
the points cu = ±\/4 + m 2 we have 

.D(cu) ~ k — 2 7 = f\/4 + m 2 7 — «(4 + ?n 2 — cu 2 ) 1 / 2 H—(4 + m 2 — cu 2 ) — z(cu =F \/4 + m 2 )y 

+ ^(4 + m 2 -cu 2 ) 3/2 + ..., 

8 

where sgn(Re \[u? + 4 — cu 2 ) = sgn(Recu). Therefore, if 7 7 ^ 0 or 7 = 0 and n 7 ^ 2 , 

(D(cu ))" 1 ~C 0 + Ci(4 + m 2 -cu 2 ) 1/2 + ..., cu ±V4 + m 2 + iO, (8.2) 
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with C 0 = (k — 2 =f i\J 4 + m 2 7) 1 and Ci = i(n — 2 =f i \/4 + m 2 7) 2 . 

If 7 = 0 and k = 2 , then in the neighborhood of the points a; = ±\/4 + m 2 , 

(^(a;))- 1 ~ i (4 + m 2 - cu 2 )“ 1/2 + J - J (4 + m 2 - cu 2 ) 1/2 + .... ( 8 . 3 ) 

2 8 

For cu —> ±m + iO , the representation (I 7 . 5 |) holds (if m/ 0 ). Hence, in the neighborhood of 
the points cu = ±m , 

D{uj) ~ k =F 27717 - i(cu 2 - m 2 ) 1/2 - z(cu =F 777)7 - i(cu 2 - m 2 ) + ^(cu 2 - m 2 ) 3/2 + ..., 

2 8 

where sgn(Re y/u 2 — m 2 ) = sgn(Recu). If m = 0 , then ( 17 . 91 ) holds. Hence, in the neighbor¬ 
hood of the point cu = 0 , 

~ 1 z 

D( cu) ~ k — zca(7 + 1 ) — -cu 2 + -cu 3 + .., cu —> 0 . 

2 8 

Hence, if (i) 7 / 0 and m 7^ 0 either (ii) 7 / 0 , m — 0 and k 0 , or (iii) 7 = 0 and 
k/ 0, then 


(D(cu))- 


_ 1 j 1 / tc + ioj (7 + 1)/k" + ..., cu —^ 0, 


if m = 0, 


Co + Ci(cu 2 — m 2 ) 1//2 + ..., cu —> ±m + iO, if m ^ 0 , 


with Cq = (k =f 77777) 1 and C\ = i(n =F 27777) 2 . If 7 = 0, m /0 and k = 0 , then 


(.D(cu)) 1 ~ i(cu 2 — m 2 ) 1//2 — ^ ^-(cu 2 — 777 2 ) 1 / 2 ..., cu —» ±777 + 2O. 


If m = k = 0 and 7 > 0 , then 


(D{u)) 


-1 


202(7 — 1) 


cu( 7 + l) 2(7+ 1) 2 8(7+ 1) 3 

Applying the inverse Fourier-Laplace formula, we have 


~h ...j (jj —y 0. 


( 8 . 4 ) 


( 8 . 5 ) 


( 8 . 6 ) 


N(t) 


+ J e-^lV(u)du, 

Im u=fi >0 


t > 0, 


with some /1 > 0. 


( 8 . 7 ) 


Note that N(t) = 0 for t < 0 . 

Theorem 8.5 Let 7, 772, k > 0 . (i) If condition C is fulfilled, then 

|IV(i)| < C(l + t)- 3/2 . 

(ii) 7/7 = 0, m 0 and k = 0 or 7 = 0 and k — 2 , then 

|IV(t)| < C(l + t)' 1/2 . 

(iii) If 7 = 0 and k > 2 , ihen 

IV(t) ~ Ci sincu 0 t + 0(t ~ 3//2 ), t —> 00 . 

(m) If rn — k — 0 , then 

N(t ) =-h 0(t _3//2 ) ; t —» cx). 

7 + 1 


( 8 . 8 ) 

( 8 . 9 ) 

( 8 . 10 ) 

( 8 . 11 ) 
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Proof (i) The proof of the item (i) is similar to the proof of Theorem 15.31 (see subsection 17.21) . 
Indeed, using Lemma [8711 we change the contour of integration in (18.71) into {a; : |a;| = R} , 


m 


i 

2tt 


e~ iujt N(u) du, 


M=* 


t > 0. 


( 8 . 12 ) 


Here R is chosen enough large such that N{u) has no poles in the region C_ D {|o;| > i?} . 
Note that if 7 = 0, then N{u) has no poles in C_ by Corollary 18.21 Denote by a 3 the poles 
of N{u) in C_ (if they exist). By Lemmas 18.11 and 18.31 there exists a 5 > 0 such that N{u) 
has no poles in the region Im u G [—5,0] . Hence, we can rewrite N(t) as 


N(t) = 


K 

-i Res a 

3 = 1 


—iujt 


N{u) I e~ iuJt N(u)du , t > 0, 

Z7T 


where £ G (0,5), the contour T mj£ surrounds of T m and belongs to the £-neighborhood of 
T m (T m;£ is oriented anticlockwise). Taking £ —* 0 , we obtain 

N(t) = j- J e~ iut ^N(u + iO) - N(u - i0)) du + o(t ~ N ) 
r m 

2 

= EE i_ f e -iut p ±^ du + o{t - N)j with any N > o 

± j =i r 

■L m 

Here P J ± (u;) := (±(u)(N(u + iO) — N(u — *0)) , where , j — 1,2, are smooth functions 

such that t(u) = 1, u G M, suppCi 1 C 0(±m) , suppC^ C C?(±\/4 + m 2 ) (Here as 
±J 

before O(a) denotes a neighborhood of a point u = a). If m = 0, we introduce Ci (-Pi) 
instead of Cf (, resp.) with suppCi C 0(0). Then (18.2ft and (18. 4 j) imply the bound (18.8ft . 
Here we use the following estimate with j = 1, 


C(u)e lu:t Va 2 — u 2 du <C(l + f) 3 ^ 2 as t —> oo, 


(8.13) 


<C(o;) is a smooth function, and ( (u) = 1 for \u — a| < 5 with some 5 > 0 (see, for example, 

[2Zj or |25]). 

(ii) Let 7 = 0, m ^ 0, k — 0. Then N(u + iO) — N(u — iO) = 2ilmN(u + iO) and N(u) has 
no poles in C_ . Introduce the circles c ±, c± = {[a; =F m\ = e} , with some £ G (0, m). We 
change the contour of the integration in (18.12)1 into T £ := U±c± U f 7 j , where 7 j , j = 1, 2, 3, 
stand for the segments of the real axis connecting the circles c± and passing in two directions, 
71 = [— y/m 2 + 4, — m — e] , 72 = \—m + £, m — e] , 73 = [m + £, \/m 2 + 4]. Using the Cauchy 
theorem and Lemma [8.11 we have 


N{t) 



c_Uc+ 


e ~ iujt N(u) du + 


3 


E 





e~ iujt lmN(u + t0)du. 
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Applying representations (j8.2j) and (18.5j) and the well-known estimate (see, for example, |25| ) 


— f e lut VuJ 2 — m 2 doj 

27r J 

|u;|=ra+l 



i cos (mt 


-tt/4 ) + 0(V 3 / 2 ), 


t —> oo, 


we obtain N(t ) = — ( 2 /( 7 r? 7 i)) 1 ' / 2 t - 1 / 2 cos(mf — 7t/4) + 0(/V 3//2 ) as f —)• oo , and the bound (18.91) 
follows. 

Similarly, if 7 = 0 and k = 2, the representations (18.31) and (18.41) imply the bound (18.91) . 
N(t ) = ( 2 /( 7 T\/m? + d)) 1 / 2 ^ -1 / 2 cos(f\/m 2 + 4 — 7t/4) + 0(t~ 3//2 ) as t —> 00 . 

(hi) If k > 2, iV(a;) has two simple poles in the points ±cuo, cuo > v ^ 2 +• Then 
calculating the rescue of e~ lut N(u)) in these points we have (I8.10p . 

(iv) Using formulas (18.61) and (j8.2[) and calculating the rescue of N(u) in the point oj = 0 , we 
obtain (18. lip . ■ 


9 Appendix C: Zero boundary condition 


Consider the following mixed initial-boundary problem on the half-line: 

{ u(x, t) = (A l — m 2 )u(x, t), igN, t G R, 

«(0,*) = 0, (9.1) 

u(x, 0 ) = uo(x), u(x, 0) = i>o(x), igN. 

In [7|, we studied the convergence to equilibrium for the harmonic crystals in the half-space in 
any dimension with zero boundary condition. However, the one-dimensional case with m = 0 
was not considered in [7j. Therefore, we outline the proof using the methods from |7j. 

We rewrite the problem (19.11) in the more general form: 

( u(x,t) = — (V{x — x')— V(x + x'))u(x',t), igN, IgR, 


x ’>0 

u( 0 , t ) = 0 , 

u(x, 0) = uq(x), u(x,0) = vo(x), x 6 N. 


(9.2) 


We assume that Wo(0) = w(0) — 0. We impose the following conditions V1-V5 on the 
interaction function V. 

VI. There exist positive constants C and (3 such that |V(:r)| < Ce~^ for 1 GZ 1 . 

V2. V is real and even, i.e., V(—x) — V{x) gR, 1 GZ 1 . 

Conditions VI and V2 imply that V(6) is a real-analytic function of 6 E T 1 . 

V3. V{6)> 0 for every 6 e T 1 . 

Let us dehne the real-valued nonnegative function, 0(0) = (U(0)) 1 ^ 2 > 0. 0(0) can be 
chosen as the real-analytic function in T 1 \ C*, where C* is a closed subset of T 1 such that 
the Lebesgue measure of C* is zero (see Lemma 2.2 in HD- 

V4. 0 /7 (0) does not vanish identically on 0 G T 1 \ C, . 

V5. (i) U(0) y 0 for all 0 G T 1 or (ii) V(0) ^0 for all 0^0, and U(0) = 0. 

Since V (0) is even, then in the second case V (0) is of order 0 2 as 0—^0. 
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Conditions V1-V5 are fulfilled, for example, in the case when V(x) has a form 


V (±1) = —1, C(x)=0 for |x| > 2, and V(Q) — 2 + m 2 with m > 0. (9.3) 


Then 0(9) = a/ 2 — 2 cos 6 + m 2 . If rn / 0 , C* = 0 , and condition V5 (i) holds. If rn — 0, 
0(6) = 2| sin(0/2)| . Then C* = {0} , and condition V5 (ii) holds. Moreover, in the case when 
V(x) is of the form (19.3ft . then the problem (19.21) becomes (19.11) . 

Write Z(x,t ) = (Z°(x, t), Z L (x, t)) = (u(x, t), u(x, t)) , Z 0 (x) = (u 0 (x),v 0 (x)). The exis¬ 
tence of dynamics is stated by the following lemma which is proved as in [7]. 


Lemma 9.1 Let atR and conditions VI and V2 hold. For any Z Q £ Tt a ,+ , there exists a 
unique solution Z(-,t) £ FL a ,+ to problem K9.2 j) . Moreover, the operator Uo(t) : Zq —> Z(-,t) 
is continuous in FL a ,+ . 


Indeed, since Gt(x) is even by condition V2, the solution of the problem (19.2ft can be 
represented as the restriction of the solution to the Cauchy problem with odd initial data on 
the half-line, 

Z i (x,t) = ^gnx-y)zi u (y), i£J, < = 0,1. 

y&L 


Here Q t 


x 


= F^ x [e A ^] with A{6) = 


(see also formulas (12.160 . (12.171 ) with 


0 1 
V(0) 0 

0(0) = (V(0)) 1/2 ) , and by definition, Z odd (x) = Z 0 (x) for x > 0 , Z odd ( 0) = 0, and Z odd (x) = 
—Z 0 (—x) for x < 0. Then the solution Z(x,t ) of the problem (19.21) is of a form 


z i (x,t) = Y,( s i i (x-y)-GHx + y))zl(y), 1 E», < = 0 , 1 . ( 9 . 4 ) 

y> i 


Moreover, the following bound holds, 

\\U 0 (t)Z 0 \\ a>+ < C( 1 + |i|) JV ||Z 0 || ai+ , (9.5) 


with some constants C = C(a), N = N(a ) < oo . 

Denote by uq a Borel probability measure on FL a ,+ giving the distribution of Zq . We 
impose conditions S1-S4 on vq . 

51. Uq has zero mean value. 

52. The mean energy density of u 0 is finite, E 0 (|Z 0 (a;)| 2 ) < C < oo. 

Finally, it is assumed that n 0 satisfies conditions S3 and S4 from Section I27Z1 

By ut, t £ M, we denote the distribution of the solution Z(x,t) = Uo(t)Zo (see Dehni- 
tion 12.61) . Then the following theorem holds. 


Theorem 9.2 Let conditions S1-S4 and VI V5 hold, and let a < —1/2 in the case V5 (i), 
and a < — 1 in the case V5 (ii). Then the following assertions are true. 

(i) The measures u t weakly converge on the space T-L a ,+ as t — y oo . Moreover, the limit 
measure is Gaussian on Ft a ,+ with correlation matrix of the form 1)2.19) - $2.20) with 

m = (vmy/ 2 . 

(ii) The correlation functions of v t converge to a limit as t —> oo, i.e., \2.1S\) holds. 
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The derivation of this theorem is based on the proof of the compactness of the measures 
family {vt,t G M} and the convergence of the characteristic functionals (cf [7]). Below we 
prove the compactness and outline the proof of the convergence for correlation functions. 

The compactness follows from the bound (19.61) below by the Prokhorov compactness the¬ 
orem [26, Lemma 11.3.1] by a method used in [ 261 Theorem XII.5.2], since the embedding 
'H a ,+ C 7-Lp,+ is compact if a > (3. 

Lemma 9.3 Let conditions SI, S2, S4 hold, and a < —1/2 if condition V5 (i) holds, and 
a < — 1 if condition V5 (ii) holds. Then the following bound holds: 

supE 0 ||[/o(t)Zo||„ + < oo. (9.6) 

t> o 

Proof. Bv Definition 12.1 1 (ii). 

Eo||2(-.()llS,+ =5L(! + |^| 2 )"(q?°(s,-) +Q”(^.-)). (9.7) 

2>0 

The representation (19.41) gives 

QlHz.z 1 ) = E 0 (z‘(z,t) ® ZHz’.t)') = (Q 0 (y,y'),Gi(y,t)®Gl,(y',t)) + , (9.8) 


where G l z (y,t) is defined in (16.10 . If condition V5 (i) holds, then the Parseval identity and 
formula (12.171) with <j>{6) = (V r (^)) 1 ^ 2 yield 

l|Gj ( .,t ) ||l = \G'-AO.t)?d0<C J (|(J(°(0)| 2 + |(J t * 1 (^)| 2 ) sin 2 (z&) dd <C 0 < oo. (9.9) 

T 1 T 1 


By conditions SI, S2 and S4, the bound (12.111) holds. Therefore, 


I Qt(z,z')\ = \{Qoi.y,y'),G l z {y,t) <g> G{,(y',t)) + \ 

< C’||Gi(-,t)|| 0l+ ||C3i(-,t)||o I+ <C' 1 <oo, 


(9.10) 


where the constant C\ does not depend on z, z' > 0 and t G M. Therefore, (19.61) follows from 
(19.71) and (19.101) . since a < —1/2. 

If condition V5 (ii) holds, then V(0) rv-' e 2 as 6 —> 0 . The estimate (19.9p with i — 1 remains 
true. Then |Q i 11 (^,P)| < C < oo by (19.101) . However, since ^ 01 (6 I ) = sin (<f>(9)t)/<f>(9), we 
have 


IGl 


< C + C l 


sin' 


fmt) 


T 1 


vie) 


sin 2 (zd)de <C + C 2 


sm 


{z9) 


sm 


0 ( 0 ) 


\0) 


dd < C + C 3 \z\, 


uniformly on ( G R. Hence, \Q®°{z,z)\ < C||G°(-,f)||o i+ < Ci + C 2 \z\ , by (19.101) . Therefore, 
for any a < — 1, 

Eol|Z(V)llh < + N 2 )“(Ci +C a N) <oo.« 

2>0 
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Now we outline the proof of item (ii) of Theorem 19.21 At first, introduce quadratic forms 
Qti'k, ’£) = (Q t (x,x'), T(:r) ® T(a/)) + , teR, $65, 


and let Q 0O ('I r , T) be the quadratic form with the kernel Qoo(x, x'). Then the assertion (12.181) 
is equivalent to the next convergence, 

lim = Qoo(d , ,'h) for any $6 5. (9-11) 

Using the operator U({t) (see formula (16.71) 1. we rewrite Q t (\ H, T) as 

Q t (tt,tt) = <2 0 (*(-, *),*(■,*)), Hx,t):=U'm= ^2^ + (z,x)^(z). (9.12) 

z&.\ 

Since Gt,+(z, 0) = Gt,+( 0,x) = 0, we put T(0) = 0, without loss of generality. Given $65, 
introduce an odd function T 0 as in Remark 12.81 Hence, 

$0M) := U'm = (9.13) 

zez 1 

Further, (see [3 p.468]) define Q*(z,z') to be equal to Q 0 (z,z') and to 0 otherwise. We split 
the function Q*(z,z') into the following three matrices (see condition S3) 

Q + (z,z') := ^q 0 (z-z'), Q~(z, z') = ^ q 0 (z - z') sign(z'), (9.14) 

Q r (z,z') := Q*(z, z') - Q + (z, z') - Q~(z, z'), zy'eZ 1 . 

Then Q t (^, VP) = (Q*(z, z 1 ), ®(z, t) <g> $(V, t)) — Y (Q a ( z , z'), $(z, t) <g> $(z',t)). Hence, the 

a=±,r 

convergence (19.111) follows from the following lemma. 

Lemma 9.4 (i) lim {Q ± (z, z'), &(z, t) <g) <F(U, £)) = {q^(z — z'), \k 0 (z) g) '& 0 (z')), with the 

£—>■00 ' 

matrices defined in \2.20 1) , (ii) lim (Q r (z, z'), <E>(z, t) (g) $(z', t)) — 0 . 

£—>•00 

This lemma can be proved using the technique of [5J Proposition 7.1]. To justify the main 
idea of the proof of Lemma [9.41 (i), we first introduce a set C 6 T 1 , 

C := C* U {e 6 T 1 \ C* : <f/(6) = 0} U {6 6 T 1 : V{6) = 0}, 

and a space S 0 := {T E S : '&<,(#) = 0 in a neighborhood of C } . Note that mes C — 0 (see (5, 
lemma 7.3]) Next, we define a norm || • ||y in the space S such that i) So is dense in S in 
this norm, while ii) the quadratic forms Q t (4/, T), 1 6 K, are equicontinuous in this norm. 
By definition, Sy is the space S endowed with the norm 

jf i (i*;(»)i 2 +(i + iv , (»)r I )i*;(e)i 2 )<«, * = 

Note that the set So is dense in , and the quadratic forms Q t (T, T), teR, are equi¬ 
continuous in Sy . Indeed, by (19.121) and (12.111) . sup |Q t (\P, \k)| < Csup ||$(-, f)||p . On the 

_ teR teR 

other hand, by the Parseval identity and (19.12)1 . (12.17ft . we have 

= [ wMfWom'de <c\\n 2 v , 

2vr Jji 
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because |^ o (0)| < C| sin#| |a;||\]/(x)| < C'(\l/)| sin#| for $65. Therefore, 

sup|Qi(T,T)| <C||T||(>, $65. 

te K 


Hence, it suffices to prove Lemma [9.41 for $ 6 So only. 

By Parseval identity and (19. 141) . 

(Q*(z, /), i(z, t) 0 i(z', i)> = 4-<»>(0), Si (9)<P»(9) 8 OTWCW). 


Note that Gt(9) = cos (j)(6)t I + sin cf>(0)t C(9). Therefore, 

1 j 


(Q + (z, z'), $(z, t ) ® $>(z', t)) = 
where Rt(9) stands for the 2x2 matrix, 


47T 


'T 1 


Rt{B)^ 0 {e)®^ 0 {e))d9, 


Rt(0) = ^[q 0 (9)+C(9)q 0 (d)C(9)*} +±cos(2mt)[U9)-C(e)q 0 (9)C(ey 
sin (2 <j>(e)t) [c{6)q 0 (e) + %(9)c(9y 


By formulas (12.20j) . one obtains 


(Q + (z, zf), *(z, t) 0 *(/, ()> = T J (g+ (9), *„(«) ® <P O (0)) d$ + ..., 

where ” ...” stands for the oscillatory integrals which contain cos(20(0))f and sin(2 4>(9))t 
with (f)(6) ^ const. The oscillatory integrals converge to zero by the Lebesgue-Riemann the¬ 
orem since all the integrands in ‘... ’ are summable by Lemma 12.41 and condition V5. The 
convergence in Lemma [9.41 (i) with sign " + " is proved. The another assertions of Lemma [9.41 
can be proved using methods of [7]. ■ 
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